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Abstract 

A large set of independent astronomical observations have provided a 
strong evidence for nonbaryonic dark matter in the Universe. One of the most 
investigated candidates is an unknown long-lived Weakly Interacting Massive 
Particle (WIMP) which was in thermal equilibrium with the primeval plasma. 
Here we investigate the WIMP abundance based on the relativistic kinetic 
treatment for gravitationally induced particle production recently proposed 
in the literature (Lima & Baranov, Phys. Rev. D 90, 043515, 2014). The new 
evolution equation is deduced and solved both numerically and also through 
a semi-analytical approach. The predictions of the WIMP observables are 
discussed and compared with the ones obtained in the standard approach. 
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1. Introduction 

In the last few years, the problem of continuous gravitational particle 
creation by a time-varying gravitational held has intensively been investi¬ 
gated in connection with the expanding accelerating Universe models [1-12]. 
The majority of the papers in this research line are somewhat based on 


Email addresses: iuribarcaiov@usp.br (I. Baranov), limajas@astro.iag.usp.br 
(J. A. S. Lima) 


Preprint submitted to Physics Letters B 


April 11, 2017 



the nonequilibrium macroscopic particle production approach discussed long 
ago by Prigogine and collaborators ini. The interest for this relativistic ir¬ 
reversible macroscopic treatment comes from the fact that it incorporates 
naturally entropy production and back reaction of the produced particles on 
the space-time geometry (see also Ref. HIES] for the associated manifestly 
covariant approach). 

Recently, we have put forward a relativistic kinetic description for the 
aforementioned non-equilibrium macroscopic treatment in the context of an 
expanding Friedmann-Robertson-Walker (FRW) geometry [16]. In partic¬ 
ular, we have argued that the covariant thermodynamic approach usually 
referred to as “adiabatic” particle production HllIS] provoked by the cos¬ 
mic time-varying gravitational held has also a kinetic counterpart. As dis¬ 
cussed there, the variation of the distribution function is associated to a 
non-collisional kinetic term of quantum-gravitational origin which is propor¬ 
tional to the ratio F/i7, where F is the gravitational particle production rate 
and H is the Hubble parameter. Therefore, for T H the process is negligi¬ 
ble and it also vanishes in the absence of gravitation regardless of the value 
of F. 

In this framework, the resulting non-equilibrium distribution function has 
the same functional form of equilibrium with the evolution laws corrected by 
the particle production process. The balance equations for entropy and par¬ 
ticle number density, as well as the temperature evolution law for “adiabatic” 
particle production were also kinetically derived both for massive and mass¬ 
less particles. 

On the other hand, there is a strong evidence for nonbaryonic dark matter 
(DM) in the Universe with a density parameter ~ 0.12 im. One of the 

most investigated candidates is a stable weakly interacting massive particle 
(WIMP) which was in equilibrium with the thermal plasma in the primeval 
Universe. In the basic scenario, one often considers a single thermally pop¬ 
ulated and coupled WIMP component whose abundance is tracked by the 
evolution equation, as phenomenologically proposed by Zeldovich [TH], and 
later kinetically derived (see, for instance, mm) 

where is the number density of the WIMP y-particles, eq hs equilibrium 
value, and < av > the thermally averaged annihilation cross section. The 
physics encoded by the above equation and even in some of its extensions 
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have already been discussed by many authors [IHl IMl HD H2] • However, as 
far as we know, neither the thermal relic abundances, nor even the extended 
form assumed by the evolution equation were discussed in the presence of 
gravitationally induced particle production. 

This paper tries to £11 this knowledge gap in two steps. Firstly, in order 
to clarify some subtleties of our previous kinetic approach [TS], we discuss 
how the particle production term modifies the above equation describing 
the evolution of thermal relics. Secondly, by assuming the popular parti¬ 
cle production rate, F^ = 3/3FT, where /3 is a dimensionless constant and 
H is the Hubble parameter, we determine how the predictions of the ob¬ 
servable quantities like the dark matter density and the thermally averaged 
annihilation cross section are modified. In the Appendix, we also show by us¬ 
ing an alternative method how the previously derived mass-shell Boltzmann 
equation for “adiabatic” particle production can be obtained (in our units 
h = ks = c = 1). 

2. Boltzmann’s Equation and Gravitational Particle Production in 

a Flat Geometry 

In principle, the simplest attempt to access gravitational particle produc¬ 
tion in a kinetic level must be expressed by some extension of the standard 
relativistic Boltzmann equation [mED] describing the evolution of the phase 
space density, /(x^,p^). A suitable generalization can be written as: 




where F((^ are the Christoffel symbols and the total collisional term, C{f) = 
Ce[f] + Ci[f], is a sum of two contributions describing elastic (e) and in¬ 
elastic (i) processes, respectively. The Vg quantity is a new source term of 
non-collisional nature encoding the particle production process due to the 
cosmic time-varying gravitational field [TB] (see also [2S])- Naturally, this 
new contribution on the r.h.s. of the above equation, has an unknown quan¬ 
tum origin whose existence is quite distinct from traditional collisional terms 
present in the standard Boltzmann treatment. 

In this work we focus our attention on the continuous particle production 
process in an expanding Universe. In this way, we consider the extended 
Boltzmann equation including particle creation in the context of a fiat FRW 
geometry: 
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( 3 ) 


ds^ = dt^ — a^{t) (dx^ + dy^ + dz^) , 
for which the non-null Cristoffel symbols are PI 

r»„ = r;j = o, = % r» =--9«. (4) 

a a 

Now, by using the mass-shell condition = m^), and the fact that 

the distribution function / is isotropic, one may show that the extended 
Boltzmann’s equation can be written as: 


E 


'dl 

dt 





(5) 


where H = a/a is the Hubble parameter, p the physical momentum and E 
is the energy of the particle. This equation was deduced in a previous work 
|16j . It generalizes the textbook mass-shell form of the standard Boltzmann’s 
equation without gravitational particle production |T9] , a process quantified 
here by the T-parameter. As discussed there, it reproduces all the thermo¬ 
dynamics results associated with the macroscopic irreversible approach to 
matter creation mms] and also incorporates the associated back reaction. 
An alternative derivation is outlined in the Appendix of the present paper. 


3. New Evolution Equation 

As discussed in the introduction, the goal of this paper is to investigate, 
based on the above kinetic theoretic approach, how the effect of gravitation¬ 
ally induced particle production changes the standard evolution equation for 
cosmological relics and the prediction of its abundance. To begin with, let 
us rewrite Eq. (|^ in the form below: 

£ol/l=C|/]. (6) 


where we have introduced the extended Liouville operator (including cre¬ 
ation) 


^c[/] 


= E 





(7) 


The parameter T^ denotes the associated gravitational production rate. Note 
that for T.^<^ H the standard Liouville operator is recovered (£c —t C). 
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dn^ 

dt 


+ 3- 


X '"x 


where 


nx = 


9x 

(2ir) 


[with dehnition (j7|)] by 


no the balance equation: 


(8) 

f fd^P, 

( 9 ) 


is the particle number density and is the number of spin degrees of freedom. 
As should also be expected, the above equation reduces to the standard one 
in the limiting case —)• 0 or, equivalently, F^ <C H. The y-particles can 
annihilate into a pair of the standard model of particle physics. For the sake 
of dehniteness, the WIMPs here are assumed to be self-conjugate Majorana 
type particles. 

Following standard lines, it is easy to see that the r.h.s. of (|^ now leads 
to an extended form of the evolution equation: 


driy^ 

dt 


+ 



< <JV > {nl-nl J), 


( 10 ) 


where is the density of the WIMP y-particles in the presence of “adi¬ 
abatic” creation (F^ ^ 0) but in absence of annihilation. In this case the 
balance equation for the number density reduces to: 


dn^^c 

dt 


T 3 '^X,C ryTT 


X"'X,c- 


( 11 ) 


Therefore, we have = n^^eq , as assumed by Zeldovich and often adopted 
in the analyzes based on Eq. Q, only if the creation rate F^ is identically 
null (this also explains our nonstandard notation). 

In principle, WIMP particles are not the unique species created by the 
time-varying gravitational held since this could violate the equivalence prin¬ 
ciple (unless the process is expressly forbidden by some unknown physical 
law). For simplicity, in what follows it will be assumed that the remaining 
particles of the thermal bath are produced by the same process. In the case 
under consideration (“adiabatic” particle production), the entropy density 
and temperature T of the thermal bath satisfy the equations below [HE]: 


- + ?,H 
s 


r 


T 

T 


a 3 ’ 


( 12 ) 
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where is the average production rate of all massless particles contributing 
to the thermal bath. 

In the standard approach (all F's = 0), the physical process is very simple. 
For T S> m^, the thermal WIMP y behaves like radiation and its equilibrium 
density is n^^eq cc T^. As the Universe expands and cools, eventually T 
and the equilibrium number density evolves like eg oc The 

concentration followed n^^eq until the density became so small that y- 
particles can no longer annihilate each other (freeze out). After that, the 
r.h.s. of (10) has no more effect on the evolution equation. Particles stop 
annihilating and lost the track from its equilibrium density, n^^eq- As we 
shall see, under “adiabatic” creation conditions, this basic qualitative picture 
is maintained just replacing n-^^eq by n^,c- However, the hnal calculated 
abundance is modihed in the presence of gravitationally induced particle 
production (see Fig. 1 and discussion below). 


4. Case of Study: F^ = F.^ = Zj3H [ 26 ] 

In order to track the evolution of the WIMP particles in a more quan¬ 
titative way, it is convenient to consider the “abundance” Y = n^/s^ where 
s is the entropy density of the thermal bath. As one may check, by intro¬ 
ducing the auxiliary variable, x = m^/T, and using the above results, the 
abundance Y satishes: 


dY 

dx 


P - P 

■^7 ^ X 


Hx 


3H J 


Y =-- 


s < av > 


Hx 


1 _Fx 

3H 


- y?) 


(13) 


This is the general relic evolution equation when WIMPs (and the remain¬ 
ing particles of the thermal bath) are being created by the time varying- 
gravitational held. As should be expected, it reduces to the standard form 


when P^ = P.y = 0 

In order to illustrate some consequences of the above equation, let us 
consider the simplest situation, namely, when the creation rate is P^ = P..y = 
3/3H [25l[26], where /I is a constant. As it appears, this kind of parametriza- 
tion is not unique. Note that for all particles the parameter (3 = T/3H 
(see Eq. (A.7)) quantihes the efficiency of the particle production process, 
and, as such, it can be a time-dependent quantity. However, by assuming 
its constancy, there is a natural upper limit, namely /? < 1 (P = 3H). In 
this case, the creation rate has exactly the value that compensates for the 
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dilution of particles due to expansion. The resulting solution of the back¬ 
ground equations is a de Sitter inflationary universe that lasts forever, that 
is, a{t) oc e^^,H = constant, p = 0 [I5]. Such spacetime resembles the 
steady state cosmology which has been ruled out by the old and present 
observation^ In addition, constant values bigger than unity also leads to 
eternal “pole” or super-inflationary expansion [27j which are also in contra¬ 
diction with the present data. In principle, one may argue that the unknown 
(5 parameter can be small since the WIMP mass can be much greater than 
the expansion energy scale at freezout. However, its precise determination 
requires an acceptable non-equilibrium theory for gravitational particle pro¬ 
duction using hnite-temperature quantum field theory in curved space-times. 
The lack of such a theory suggests naturally a phenomenological approach. 

In what follows, it will be assumed that is constant and smaller than 
unity. For this toy model, we see from Eqs. ( [IT| ) and (12) that the particle 
density and entropy in the absence of annihilations are: 


n^^c oc a 


-3(l-/3) 


s OC a 


-3(l-/3) 


(14) 


These results suggest that Y must evolve to a constant value at very late 
stages. We now demonstrate this result by integrating directly the evolution 
equation. 

A considerable part of the process occurred during the radiation-dominated 
era. In addition, if the creation happens “adiabatically”, that is, in such a 
way that the specihc entropy (per particle) remains constant, it means that 


the temperature evolves obeying (12), but the entropy and energy densi¬ 
ties take the equilibrium forms [26l l25] : s(T) = p{T) = 

where and g^, counts for different relativistic degrees of freedom (d.o.f.) 
Here we neglect any temperature variation of the relativistic d.o.f. by tak¬ 
ing g = g*,s = S'* = 90 izu. As one may check, using these results and 
dehnitions, the evolution equation becomes: 


dY 

dx 


A 


(y" - Y'l) , 


(15) 


(1 — jd) x'^ 

where A = 0.264m^mp;y^ < av >. For typical values ~ 100 GeH and 


^Note that such a solution is not supported by the standard vacuum energy density. 
For instance, we may have a de Sitter solution with p = p/3. 
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Figure 1: Solutions of the new evolution equation for different values of (3. The product of 
the mass by < cru > was fixed such that the parameter A = 10^^ (see its definition below 
Eq. (15l). For a given value of m^, the effect of the (3 parameter is fully equivalent to 
increasing the cross section. 


< av >~ 3 X 10 ® GeV one finds A ~ 10^^. Note that for /3 = 0 the 
standard result is recovered |2I- 

In Figure we display the predictions of the new evolution equations 
for distinct (3 values. The freeze-out of the WIMPs is delayed thereby lead¬ 
ing to a lower /^-dependent abundance in comparison with the standard case 
(/? = 0). As remarked before, the overall evolution picture is maintained in 
the presence of gravitationally induced particle production. This happens 
because for an “adiabatic” particle production process the temperature evo¬ 
lution law is modihed, but the equilibrium form of the phase space density 
is preserved [2Sl [IS] • 

The contrast between the standard and the extended approach as dis- 

Y-Y j. 

cussed here, can also be quantihed by percentage difference, AY = y , 

st 

where Ygt is given by the standard equation without creation (/9 = 0). 

In Figure 2 we plot this percentage as a function of x for different values 
of (3. As one may see, the relic abundance can be changed by large factors 
(10% to 30%) which depend on the values assumed by the (3 parameter. 

As usually, the associated cosmological parameters can also be extracted 
even based on a semi-analytic solution [201 El]- At early times, x ^ Xf, where 
X/ is the X at freeze-out, the exact solution of the evolution equation follows 
the “equilibrium solution” closely. Introducing a new variable, A = Y — Yc, 
the new evolution equation takes the form 




A A 


(2Fc + A), 


dA 

dx 


dx (1 — (3) 


(16) 












Y — Y j- 

Figure 2: Percentage ratio AY = y for different values of /3. As in Fig. 1, for all 

st 

curves we have fixed the parameter A = 10^^ (see text). 


which has the following solution for x !:§> 1 


(l-/3)x2 

2A 


(17) 


where it was used that ^ ~ —Y, A 1 and dA/dx 1 in this regime. 

Freeze-out occurs when A = sYc, where £ is of order unity. Putting this 
condition into the above result, one obtains an equation for Xf 




In 


0.156em^mpi < av > 
( 1 -/ 3 )^/^ 


(18) 


using 5 = \/2 — 1 [21] and the assumed typical WIMP values, it is readily 
checked that the freeze-out occurs at x/ ~ 21 — 22, depending on the value 
oi(3. 

At very late times (x !:§> x/), the efficiency of relic production by the 
thermal bath becomes negligible, this means that the term < av > can 
be ignored and the equation to be solved then is 


dY _ _ A 
dx (1 — /9) x^ ’ 


(19) 


which has the solution 

F(x —)■ cxd) = -^x/. 

A 


( 20 ) 


Remembering that = m^SoY{oo) and using sq = 2900cm“^, 

one can verify the effects of (3 on the function of < cm >. More 

precisely. 


9 

















6.x10“'“ l.xlO"® 1.4»10"®1.8»10“^ 


<crv>[Ge\r^] 


Figure 3: Relic density as function of < cru >. The shaded region represents the 

constraints = 0.1196 ± 0.0031 (Icrc.L) as provided by the Planck experiment m- 


Q h'^ = = {1 - (3)8.5 X 10-^^ - ^ -. (21) 

The expressions (18), (20) and (21) can be compared to Equations (10), 
(11) and (13) in [21]. In absence of particle production (in this case, (3 = 0) 
and using their expansion for < av >, both results are equivalent. 

In Figure 3, by assuming ~ lOOGeE, we display the behavior of 
as a function of the thermal averaged cross-section for different values of (3 as 
depicted in the hgure. The shaded region delimits the upper and lower limits 
as given by the latest constraint, = 0.1196 ± 0.0031 (lac.Z.), provided 
by the Planck experiment im. The solutions have been plotted both for the 
numerical (solid lines) and the semi-analytical treatment (dashed lines). 

In Figure 4, by adopting again the above Planck constraint on 
we show how the thermal averaged cross section < av > depends on the 
(3 parameter. As a general rule, for a given we see that the value of 

< av > decreases for increasing values of the (3 parameter. 

5. Conclusions 

In this paper we have deduced a new evolution equation for thermal relics 
in the context of a relativistic kinetic approach for particle production re¬ 
cently proposed [16]. By assuming a flat FRW cosmology, the equation was 
numerically solved for a creation rate P^ = ?>(3H, where /3 is a positive con¬ 
stant parameter. The results can be summarized by the following statements: 

(i) For all physically meaningful values of (3, the basic qualitative picture is 
maintained and the hnal abundance saturates when the annihilation becomes 
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Figure 4: Annihilation cross section < crn > as function of the /3 parameter. Like in the 
previous Fig. 3, the lines also delimit the Itr upper and lower bounds of from Planck 
collaboration. 


negligible. The final calculated abundance is modified in the presence of 
gravitationally induced particle production (see Fig. 1). 


(ii) For the adopted creation rate (F^ = 3(3H), the new evolution equation 
differs only by a pre-factor of 1/(1 — /3) as compared to the standard case 
(see Eq. ([I^). 


(hi) The percentage difference relic abundance in the presence of creation 
can be changed by large factors (10% to 30%) which depend on the values 
assumed by the [3 parameter (see Fig.2). 


(iv) The present day cold dark matter density and the thermally averaged 
annihilation cross section are also strongly dependent on the (3 parameter 
(see Figs. 3 and 4). For all studied values of /3, the thermally averaged 
cross section should be smaller than the usual case (when /3 = 0) in order to 
accommodate the latest Planck results. 


Finally, we stress that the new evolution equation and the toy model dis¬ 
cussed here open the possibility of more realistic analyzes either by assuming, 
for instance, a more plausible matter creation rate, the variation of the effec¬ 
tive relativistic d.o.f. or even the possibility of a WIMP decay process after 
the “freeze-out”. Some discussions along these lines will be presented in a 
forthcoming communication. 


Appendix A. On the Kinetic Equation with particle production 

In this Appendix, we provide an alternative derivation of the mass-shell 
Boltzmann equation for a fiat FRW geometry in the presence of gravitation- 
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ally induced particle production [see Eq. ([^]. To begin with let us consider 
the extended covariant Boltzmann’s equation in the form [16] 

where for simplicity we have neglected the collisional contributions C{f). The 
main goal here is to model the new quantity of non-collisional nature, Vg, 
associated with the particle production process. In principle, the intensity 
of the process depends on the scalar ratio T/G, where 0 = = 3H is 

the scalar of expansion (macroscopic time scale). Since Vg must take into 
account gravity, the simplest ansatz is [T6] : 


•p =_ 

'9 0 a/sP 


p^L 


91 


(A.2) 


where Lg is an unknown 4-vector that can be determined by considering the 
homogeneous conditions of the FRW metric where the covariant irreversible 
thermodynamic approach was originally formulated [m dsi US] • Elementary 
considerations implies that Lg can be written as a linear combination 


^9 ^Qpfi ^^91 (^-3) 

where A and B are arbitrary constants and = 0. By using this con¬ 

straint and evaluating the result in the comoving frame one Ends: 


Lg = A 


df _ Pm 

dpf^ Po dp^ 


and inserting this result into (|A.2|), the Boltzmann equation (|A.1|) 

.9^ 




dx9 




fjP P 




dp- 


-9 


0 


df _ Pg df 

dp9 Po dp^ 


(A.4) 

becomes 

(A.5) 


How to fix the constant A? By using the kinetic dehnition of the particle 
flux a simple algebra shows that A = 1/2 in order to recover the macroscopic 
approach (for details see section IV in [IS]). Now, by adding the collisional 
term, the above equation in the mass-shell becomes 


4 = 


(A.6) 
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Finally, by introducing the physical momentum (p) which is related to the 
comoving momentum (p) by p = ap, after a simple change of variables the 
Boltzmann equation takes the form: 

The above equation should be compared with Eq. (§ in section II or Eq. (47) 
appearing in our previous article |I6]. Note that in Eq. (|^ we have adopted 
physical quantities but the bars were dropped as usual in the literature. 
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